Role of parton shadowing in the comparison of p-A and A-A results 
on J/^ suppression at energies available at the CERN Super Proton Synchrotron 
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The observation of an anomalous J/i/; suppression in nucleus- nucleus collisions is one of the most 
important results of the SPS heavy-ion program. An essential ingredient in this result is the determi- 
nation, obtained by studying p-A collisions, of effects not related with the formation of a deconfined 
medium. These effects are extrapolated to A-A collisions, determining a reference J/t/j yield which 
is then compared with the measurements. In this article we investigate the role of parton shadowing 
on the determination of such a reference, and we calculate its effect for In-In and Pb-Pb collisions 
as a function of rapidity and centrality. 

PACS numbers: 25.75.-q, 14.40.-n 



I. INTRODUCTION 

The suppression of charmonium states was proposed a 
long time ago as a signature of the production of a de- 
confined state in nucleus-nucleus collisions [1]. However, 
it was very soon realized that nuclear effects not related 
to deconfinement or, more generally, to the production 
of a hot medium may influence the observed charmonia 
suppression 0, |3| • Such effects were investigated through 
the study of charmonia production in p-A collisions. In 
these reactions, the produced cc pair may interact with 
the cold nuclear medium of the target nucleus, hindering 
the formation of a bound state. 

Several p-A data samples exist today for J/?/; produc- 
tion at fixed target energies, in particular from NA50 
at the SPS from E866 at FNAL [Z] and from the 

HERA-B Collaboration at HERA Nuclear effects are 
usually parametrized by comparing the yields for vari- 
ous nuclear targets in a certain kinematical region, and 
then fitting their A dependence in terms of the simple 
power law A°' . Alternatively, the data are analyzed in 
the framework of the Glauber model @, and their A 
dependence expressed as an effective "absorption" cross- 
section (Tj^^. The interpretation of J/V' results in p-A II. SHADOWING EFFECTS: p-A COLLISIONS 

collisions is the object of a rather strong theory effort. 
Nowadays, it is clear that not only the nuclear disso- 
ciation of the cc pair plays a role, but also effects like 
shadowing, initial- and final-state parton energy loss, and 
possibly the intrinsic charm component of the projectile 
should be taken into account in a realistic description of 
the results [l3|. 

When studying i/ip suppression in A-A collisions, a 
precise knowledge of nuclear effects is an essential requi- 
site to disentangle genuine hot-medium effects. In the ap- 
proach commonly used up to now [sj , the effective quan- 
tity cTjy^ is obtained analyzing p-A data taken in the 
same kinematic domain of A-A collisions under study. 
Then, it is assumed that in both p-A and A-A collisions 



scale with L, the mean thickness of nuclear matter seen 
by the cc pair in its way through the projectile and tar- 
get nuclei. In this way it is possible to determine the 
expected JM yield for nuclear collisions as a function of 
centrality With this approach, a significant anoma- 
lous suppression (i.e., a suppression which goes beyond 
the estimated contribution from nuclear effects) has been 
detected at SPS energies [llllTl|. 

In this article, we go a step further by taking explicitly 
into account parton shadowing in the determination of 
the cold nuclear matter effects in A-A collisions. Shad- 
owing of partons in nuclei is a depletion of their pop- 
ulation at small momentum fraction of the nucleon, x, 
compared to that in a free nucleon, with a corresponding 
enhancement at moderate x (anti-shadowing) . Contrar- 
ily to final state dissociation of the cc pair, shadowing is 
not expected to scale with L, because in p-A collisions 
only partons in the target are affected by shadowing, 
whereas in A-A the projectile is also involved. There- 
fore, a different approach to the evaluation of nuclear 
effects is required. 



To give an estimate of shadowing effects on 3 /if) pro- 
duction, we use the Color Evaporation Model (CEM) at 
leading order (LO) [l^- In this approach the charmo- 
nium production cross section for p-A collisions is ob- 
tained by integrating the free cc cross section from en- 
ergy threshold to the open charm threshold. In absence of 
final-state interactions of the produced J/f/' one has [l3| 



_W ^ 2F / mdm 

A dXp J2me 



gp^(xi,X2,TO^) 



(1) 



nuclear effects, parametrized through the quantity cTj 



where F is the fraction of cc pairs which gives a J/'0 
in the final state and HpA is given by 
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HpA{xi,X2,m'^) = f^{xi,m^)f^{x2,m^)(Tgg{m^)+ 
+ f^{xi,m')f^{x2,m^)a,g{m^) 

q—u.d,s 



(2) 



is the sum of two terms, corresponding to the 
elementary cc production by gluon fusion and qq anni- 
hilation, convoluted with the parton densities ff{x,Q^) 
and /f (x, Q^) in the target and projectile nucleons, eval- 
uated at Q = m. xi {X2) is the fraction of the nucleon 
momentum carried by the parton of the projectile (tar- 
get) which interacts. We have used in our calculation 
the GRV98LO [l^ set of parton distribution functions 
(PDF). The PDF for nucleons inside nuclei are modified 
with respect to the free ones according to the following 
expression: 



ft{x,Q-')=Rt{x,Q-')f!'{x,Q-') 



(3) 



Various parametrizations of the PDF modifications 
from nuclear effects exist. We have used the EKS98 [l^ 
and EPS08 [l^ sets, which are available for all mass num- 
bers and have been implemented in the frame of the com- 
monly used LHAPDF interface [l^ . The calculation has 
been performed for several A values, corresponding to 
Be, Al, Cu, In, W, Pb and U target nuclei. For each tar- 
get the ratio between the J/ip production cross section 
(per nucleon-nucleon collision) in p-A and p-p, which we 
refer to as shadowing factor in the following, is given by 



(xf) 



1 dcT^j/^/'^^F 
A daY/^/dxp 



(4) 



In Fig. [T] we present the result of the calculation of 
the shadowing factors for p-A collisions at 158 GeV, the 
energy used in the A- A data taking at the SPS. The 
plots refer to midrapidity, and we separately show the 
shadowing factors for the gg and qq fraction of the pro- 
duction cross section. We plot the results as a function 
of the L variable fl9||, computed for each nucleus using 
the Glauber model with realistic nuclear density distri- 
butions [20l |. 

It can be noted that both EKS98 and EPS08 sets give 
shadowing factors larger than 1, and that such an anti- 
shadowing is more pronounced when using EPS08. The 
qq fraction of the cross section exhibits a small shadow- 
ing, which is more than counterbalanced by the dominant 
gg production process. We have also verified, by choosing 
another set of PDFs (MRST2001LO [U) that our results 
do not depend, within at most 2%, on the specific choice 
of the PDF set. The flattening of S'^JJ' for heavy nuclei 
when using EPS08 is because of the fact that a constant 
shadowing is assumed beyond A=208 [l^. In Fig. [2] we 
present the results of a similar calculation at slightly for- 
ward {y = 0.5, y = 1) and backward {y = —0.5, y = — 1) 
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FIG. 1: Shadowing factors for pA collisions at 158 GeV, at 
midrapidity. Open circles refer to the gg fraction of the cross 
section, open squares to qq, solid circles to the total cross 
section. The top plot has been obtained with the EKS98 set, 
the bottom plot with EPS08. 



rapidity. We also plot the y — result for comparison. 
We recall that the region < y < 1 is the one where J/ip 
results are available at SPS energies for A-A collisions. 
The shadowing factors are found to depend on y, as ex- 
pected, since X2 is directly related to this quantity. We 
recall that X2 — (toj/^,/^/s) exp(— y) for J/ip production 
in the GEM at LO. In Fig. [3] we plot the EKS98 and 
EPS08 parametrizations of nuclear effects on PDFs for 



the Pb nucleus, relative to Q 



The arrows indi- 



cate the X2 values corresponding to the rapidities where 
we have performed the shadowing factor calculation of 
Fig. [21 It can be clearly seen that from y = —1 towards 
y = 1 one goes from a shadowing to an anti-shadowing 
regime, reaching a maximum of the anti-shadowing effect 
at X2 --0.10 for the EKS98 set (~0.13 for EPS08), which 
then decreases going towards smaller X2 (corresponding 
to more forward rapidities). Such an evolution can clearly 
be identified looking, at fixed L, to the shadowing factors 
shown in Fig. [21 which are <1 at y = — 1, then increase 
and finally decrease at forward rapidity. 

We note that the shadowing factors calculated for 
y = —0.5,-1 are identical, for symmetry reasons, to 
those for A-p collisions at y = 0.5,1. Results on A-p 
collisions are interesting in order to study effects from 
PDF modifications in the projectile, which will be im- 
portant when studying A-A interactions. We see, for 
example, that at y = 1 the shadowing effects in A-p are 
very different with respect to p-A. 
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FIG. 2: Shadowing factors for J/^/; production in p-A colli- 
sions at 158 GeV, using EKS98 (top plot) and EPS08 (bottom 
plot). Solid and open triangles refer to j/ = 1 and y = —1 re- 
spectively, solid and open squares to y = 0.5 and y = —0.5, 
circles to j/ = 0. 



FIG. 3: The EKS98 (continuous line) and EPS08 (dashed 
line) parametrizations of nuclear modifications to the PDFs 
for the Pb nucleus, calculated at Q = mj/^. The arrows 
correspond to the rapidity values where the calculation of 
Fig. [2] has been performed. 



III. SHADOWING EFFECTS: A-A COLLISIONS 

When moving to A-A collisions, the cross section and 
the shadowing factors, integrated over centrality, can be 
calculated with an expression similar to Eq. [U by replac- 
ing HpA(xi,X2,rn?) with if^^(a;i, X2, m^), taking now 
into account that shadowing affects both projectile and 
target nuclei. The problem becomes more complicate if 
shadowing factors have to be calculated for various cen- 
trality intervals. Clearly, for various geometries of the 
collision, either the halo or the core of the nuclei will 
be mainly involved, and the shadowing effects will be 
more important in the core than in the halo. Various 
parametrizations of the local shadowing inside the nu- 
cleus have been proposed. We have used two of them [2^ : 



Rt^{xM\r,z) = l + Nf{Rt{x,Q') (5) 



and 



/ dzpA{0,z) 

In the first one, shadowing in a certain location (r, z) 
inside the nucleus is proportional to the local nuclear 
density PAi'f'^ z), while in the second it is proportional to 
the length L of nuclear matter crossed by the parton on 
its way through the nucleus. The normalization is 



fixed to ensure that / drdzRf_p{x,Q^ ,f, z) = Rf{x,Q^) 
(and similarly for Nf). 

The study of the shadowing factors in A-A collisions 
as a function of centrality has been performed for In-In 
and Pb-Pb at 158 GeV/nucleon (see also [ISl for a pre- 
vious investigation of the infiuence of shadowing on the 
centrality dependence of J/ip and Drell-Yan cross sec- 
tions). A large number of events has been generated for 
centrality values in the interval 0< b <12 fm for In-In 
and 0< b <16 fm for Pb-Pb, in steps of 2 fm, with a 
Glauber Monte-Carlo approach. We have used o'f^g;=30 
mb and the measured nuclear density distributions for In 
and Pb 20] . For every N-N collision in each A-A interac- 
tion, we calculate pA{r, z) and L(f) for the two colliding 
nucleons, and the product of the two corresponding shad- 
owing factors, according to Eq. [5] and [6l By averaging 
the shadowing factors over all the N-N collisions in each 
A-A interaction, we get the centrality dependence of the 
shadowing factors. In Fig. |4] we show, as a function of 
i, the calculated shadowing factors for In-In and Pb-Pb 
collisions at y—Q.b. For the SPS data at 158 GeV, this is 
the rapidity where the acceptance reaches its maximum. 
The symbols, connected by the continuous lines, have 
been calculated using Eq. [5] for the local dependence of 
shadowing inside the nucleus, whereas the dashed line has 
been obtained using the parameterization of Eq. [6] We 
notice that the two parametrizations give similar results, 
their difference not exceeding 2-3%. In absence of other 
nuclear effects, this result implies a ^-^10% anti-shadowing 
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effect for central nucleus- nucleus collisions. When using 
EPS08, the effect increases up to ~25%. 
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FIG. 4: Shadowing factors for In-In (top) and Pb-Pb collisions 
at 158 GeV, at j/=0.5, obtained with EKS98. Solid circles re- 
fer to the EKS98 parametrization, open circles to EPS08. The 
symbols connected by a continuous line have been calculated 
using Eq. [5] whereas the dashed lines have been obtained us- 
ing the parametrization of Eq. [5] for the local dependence of 
shadowing. 



IV. COMPARISON OF COLD NUCLEAR 
MATTER EFFECTS IN p-A AND A-A 
COLLISIONS 

We now want to discuss the extrapolation of cold nu- 
clear matter effects measured in p-A to A-A collisions, 
taking into account shadowing effects together with final 
state absorption of the created cc pair. To do that, we 
model final state absorption effects using the simple law 
exp(— pcTj^^L), which has been shown to reproduce, at 
first order, data on i/il) production at fixed target ener- 
gies. In this approach, the i/ip production cross section 
per N-N collision in p-A is given by 



A 



pA 



NN 



S'J^ X eM-p^f/%L) (7) 



We then fit the p-A cross-section per N-N collision, ob- 
tained through Eq.[7]with the simple exp(— pcTjy^L) law, 
i.e. neglecting the existence of shadowing. This is the 
procedure usually followed at SPS energies to extract the 
so-called "normal nuclear absorption" . Clearly this ap- 
proach, when shadowing factors are not negligible, leads 



to two main consequences. On one hand the obtained 
(TjJ^ values do not represent anymore the size of final- 
state absorption, but an effective quantity which, due to 
the presence of shadowing [25l [26j , may be quite different 
from CTj^^. On the other hand, because shadowing effects 
do not necessarily have the same L-dependence of final 
state absorption, exp(— pCj^^L) may not give anymore a 
reasonable fit of the p-A data. However, given the size of 
the shadowing corrections, their deviation from the ex- 
ponential behavior is difficult to observe in the existing 
data samples. 

The i/i}} cross section per N-N collision in A-A, as a 
function of centrality, is then obtained using Eq. [71 re- 



placing S^'^ with SjI^ and \/A with l/iVcon, where the 
shadowing factor and the number of N-N collisions are 
calculated for the various centrality bins. We can now 
compare the extrapolation of the exponential fit of the 
p-A calculation with what has been obtained for A-A col- 
lisions. When doing that, we may expect significant de- 
viations from such an extrapolation. First of all, in A-A 
shadowing affects not only the target nucleus but also 
the projectile, leading to an extra-effect that is clearly 
not present in the naive extrapolation of the exponential 
fit of p-A. Furthermore, these deviations may heavily de- 
pend on rapidity, because when moving away from y=0 
the a;-region probed in the projectile becomes different 
from that of the target. 

To illustrate the procedure, we show in Fig. [51 for 
158 GeV energy and j/=0.5, the expected behavior for 
p-A collisions, including shadowing and having assumed 
cr y^=4: mb. We also show the result of an exponential 

fit to p-A data, which gives tTjy^=3.0 mb, a value smaller 
than (Tj/^ , because of the presence of an anti-shadowing 
effect. In the same plot we compare this exponential fit 
with the calculated cross section values per N-N colli- 
sion for In-In and Pb-Pb, as a function of centrality. We 
clearly see that A-A cross sections deviate from the ex- 
trapolation of p-A results. In particular, at fixed L, the 
J/V' cross section per N-N collision is systematically lower 
in A-A with respect to p-A. However, we note that the 
relative behavior of A-A results with respect to p-A, at a 
certain L, cannot be easily deduced by a simple inspec- 
tion of the p-A (and A-p) shadowing factors. In fact, 
from the geometry of the interaction, it can be shown 
that the same L for p-A (A-p) and A-A corresponds to 
very different average nuclear densities probed in the col- 
lision and therefore to a different average strength of the 
shadowing effects. (To give a numerical example, the L 
value corresponding to p-Pb collisions is obtained with 
Pb-Pb collisions at = 12 fm. The average nuclear den- 
sities probed are 0.76po and 0.48po respectively, where 
po is the core nuclear density). Therefore, the A-A shad- 
owing factors cannot be obtained as a simple product of 
p-A and A-p shadowing at the same L. 

A deviation of A-A results from the p-A extrapola- 
tions was indeed found at SPS ener gies , and it was called 
"anomalous J/^jj suppression" [iTI. K^. Usually this ef- 



5 



0.5 



III 
8 10 
L(fm) 



c 1.1 

o 

w 

$1.05 

a. 

a. 

ID 1 

= 
O 

■(5 0.95 

E 
o 

< 0.9 



0.85 
0.8 

0.75 
0.7, 



: ■ • ■ • 
■ 




y=0 
■ ■ ■ 


■ 

I ■ • 

■ 

■ 

— • 

■ • 


■ 

• •• 

■ 


y=0.5 

■ ■ ■ 

y=1 

■ ■ ■ 


-i 1 , , 1 , , 1 , 1 , 







50 100 150 200 250 300 350 400 450 

''''part 



FIG. 5: Comparison of p-A and A-A cross sections for J ftp 
at 158 GeV, for y=0.5, including shadowing and final state 
absorption, in arbitrary units. The p-A cross sections per 
N-N collisions are shown as triangles. The hne represents an 
exponential fit to p-A data. In-In and Pb-Pb cross sections 
are shown as circles and squares, respectively. 



FIG. 6: Ratio between the A-A cross section per N-N collision 
and the extrapolation of p-A results, for three rapidity values, 
using the EKS98 parametrization. The results for In-In and 
Pb-Pb are shown as circles and squares, respectively. 



feet was connected to hot nuclear matter effects, includ- 
ing the production of a deconfined state [27] . The result 
of Fig. [5] shows that at least a fraction of this effect can 
be attributed to having neglected the influence of shad- 
owing in the determination of the "nuclear absorption" 
reference. In Fig. [6] we present the ratio between the 
J /ip cross sections for In-In and Pb-Pb and the exponen- 
tial extrapolation of p-A results, for three rapidity values 
{y=0, 0.5 and 1), using EKS98 and Eq. [5] for the local 
dependence of shadowing. This kinematical range cor- 
responds to the region where 3/^ production has been 
studied by the NA50/NA60 experiments in nuclear col- 
lisions. The values of these ratios do not depend on the 
specific o-j^^ value used in the calculation, because, at 

fixed L, the factors exp(— pUj^^L) cancel out in the ratio 
of the cross sections between A-A and p-A. 

The values plotted in Fig. ^ show that a simple ex- 
trapolation of p-A results, obtained through a fit of the 
A-dependence that parametrizes all nuclear effects under 
a single parameter CTjy^, is not in agreement with A-A 
3/ip cross sections as a function of centrality. This is 
basically because of the presence in A-A of shadowing 
effects in the projectile that are of course not present in 
p-A. Although the effect remains rather small at midra- 
pidity, it increases fast when moving away from y—0, 
reaching a discrepancy of ^^20% between the extrapo- 
lation of p-A data and the central A-A results for y=l. 
This effect must clearly be taken into account when look- 
ing for hot nuclear matter effects in A-A data. In partic- 



ular, the anomalous suppression values observed at the 
SPS must be rescaled by the values shown in Fig. [6] at 
the corresponding rapidities, resulting in a ~10% aver- 
age reduction of this effect in the interval 0< y <1 where 
data were taken. When using the EPS08 parametriza- 
tion, similar values to those of Fig.[S]were obtained, with 
relative discrepancies not larger than ~2%. However, it 
is well known that the uncertainties on the modification 
of gluon PDFs are quite large, and very recent analysis 
(EPS09) arc now starting to systematically address this 
issue 30]. The LO set of EPS09 gluon PDFs modifica- 
tion has an average value, in the a;-region corresponding 
to SPS data, quite similar to the EKS98 one, with an 
error of the order of ±15%. By injecting such an un- 
certainty in our calculation, it turns out that the values 
shown in Fig. |6] vary by about 5%. Finally, other sets of 
nuclear modifications to gluon PDFs exist in the litera- 
ture that exhibit either no (or little) anti-shadowing in 
our X- region (nDS/nDSg [28*]) or an anti-shadowing very 
strongly increasing with x (HKN [1^). The use of such 
sets in our analysis gives almost no difference between 
the extrapolation of p-A results and A-A for nDS/nDSg 
(ratios ~1 in Fig. [B]), or higher values, increasing with y, 
for A-A with respect to p-A for HKN (i.e. ratios larger 
than 1 in Fig. IH]). However, it was pointed out pol] that 
such analyses might be less constrained in the a;-region 
under study, since they do not make use of data from 
high-pT TT production in d-Au collisions from the BNL 
Relativistic Heavy-Ion Colhder (RHIC) [H [H, which 
are relevant for the determination of the large- a; region 
gluon contribution. 
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V. CONCLUSIONS 

We have investigated in this article the role of shadow- 
ing for i /ij} production in p-A and A- A collisions at SPS 
energies. In particular, we have shown that an extrapola- 
tion of cold nuclear matter effects measured in p-A which 
does not take explicitly into account shadowing, cannot 
reproduce in a correct way such effects for A-A. In the 
frame of an LO Color Evaporation Model calculation, 
performed using the EKS98 and EPS08 parametriza- 
tions, we have shown that neglecting shadowing, the p-A 
extrapolation is biased by ^10% at y=:0.5 for central 
In-In and Pb-Pb collisions. Such a bias must be taken 
into account in analysis that aim at determining effects 



from hot nuclear matter m i /ip production in A-A colli- 
sions at SPS energies (the so-called "anomalous suppres- 
sion"). In particular it may be quantitatively important 
for lighter systems such as In-In, where the deviations 
from p-A extrapolations are relatively small. 
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